Polarizabilities, Stark shifts, and lifetimes of In atom 
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We evaluate the polarizabilities of the 5pi/2, 6s, 6pi/2, and 6^3/2 states of In using two different 
high-precision relativistic methods: linearized coupled-cluster approach where single, double and 
partial triple excitations of the Dirac-Fock wave function are included to all orders of perturbation 
theory and an approach that combines the configuration interaction and the coupled-cluster method. 
Extensive comparison of the accuracy of these methods is carried out. The uncertainties of all 
recommended values are evaluated. Our result for the 6s — 5pi/2 Stark shift is in excellent agreement 
with the recent measurement [Ranjit et al, arXiv: 1302.0821 vl]. Combining our calculation with 
this precision measurement allows us to infer the values of the 6pi/2 and 6^3/2 lifetimes in In with 
0.8% accuracy. Our predictions for the 6^3/2 scalar and tensor polarizabilities may be combined 
with the future measurement of the 6s — 6P3/2 Stark shift to accurately determine the lifetimes of 
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I. INTRODUCTION 

An indium atom represents an excellent system to com- 
pare the accuracy of different high-precision theoretical 
methods since it may be considered as both the monova- 
lent system (assuming a closed 5s 2 core) and a trivalent 
system. In this work, we use both linearized coupled- 
cluster approach where single, double and partial triple 
excitations (LCCSDpT) of the Dirac-Fock wave func- 
tion are included to all orders and a method combining 
the configuration interaction (CI) and the coupled-cluster 
method (CI+all-order) . 

The monovalent LCCSDpT method has been applied 
to a large number of neutral and ionized monovalent sys- 
tems and yielded very accurate predictions for a num- 
ber of atomic properties (see review [l[ and references 
therein). It has been used for a variety of applications 
ranging from the study of fundamental symmetries 0, H[ 
to study of the degenerate quantum gases and quantum 
information Q. This method has been tested against 
all high-precision alkali and monovalent ion experimental 
values which allowed to establish a systematic procedure 
to evaluate its uncertainties even when no experimental 
data exist 

The CI+all-order method was recently developed for 
the treatment of the more complicated systems m. It 
has been tested on a variety of divalent systems [7Hllj 
and applied to Tl [12|. However, there are far less ex- 
perimental benchmark experiments for the properties of 
the divalent and trivalent systems in comparison with the 
alkalis. A recent high-precision (0.27%) measurement of 
the 6s — 5pi/2 Stark shift in In [13j | provides an excellent 
opportunity to compare the accuracy of the CI+all-order 



and LCCSDpT approaches, and test the procedure for 
the evaluation of the uncertainties of theoretical data. 
In turn, the combination of the Stark shift measurement 
[13j and our calculations allows us to infer the values of 
the 6pi/2 and 6p 3 / 2 lifetimes in In with 0.8% accuracy. 
It is extremely difficult to achieve such high accuracy via 
direct lifetime measurements. Moreover, we predict the 
values of the 6p 3 / 2 scalar and tensor polarizabilities which 
can be soon tested against the future measurement of the 
6s — 6^3/2 Stark shift [l3| and will allow to determine the 
5dj lifetimes. 

Precise knowledge of In properties is of interest to the 
study of the fundamental symmetrises, including par- 
ity violation and search for the permanent electric-dipole 
moment, since it is similar to Tl (also group IIIB). Un- 
derstanding of the theoretical and experimental uncer- 
tainties in In can be used in Tl studies. Both EDM [l4| 
and parity- violation studies [HI, [r| had been carried our 
in Tl. The theory accuracy in Tl is still below the ex- 
perimental accuracy, hindering further parity violation 
studies with this system. A recent controversy regarding 
the calculated values of Tl EDM enhancement factor is 
discussed in detail in (l2l |. 

In-like ions are excellent candidates for a search for the 
variation of the fine-structure constant a. Despite very 
large ionization energies, certain ions have transitions 
that lie in the optical range due to level crossing and are 
very sensitive to a- variation [17]. In-like ions are partic- 
ulary well suite for the experimental search for such tran- 
sitions 18]. In fact, In-like isoelectronic sequence has by 
far the largest number of ions with long-lived metastable 
states with transition frequencies between 170 and 3000 
nm, high sensitivity to cv-variation, and stable isotopes 
[l8| . One of the main obstacles for the experimental work 
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in this direction is the lack of any experimental data for 
these systems and difficulty of accurate theoretical pre- 
dictions. Testing the CI+all-order method on neutral In 
provides important information of the accuracy of this 
approach for further studies of the In-like ions. 



II. METHOD 

We use two different relativistic high-precision ap- 
proaches for all calculations in this work. Such com- 
parison of these two methods for the same system have 
never been done before. The first method (LCCSDpT) 
has been extremely successfully in predicting properties 
of alkali- metal atoms and other monovalent ions [l|. It 
can also be applied to In by treating 5s 2 shell as a part of 
the [ls 2 2s 2 2p 6 3s 2 3p 6 3(i 10 4s 2 4p 6 4ci 10 5s 2 ] core. The disad- 
vantage of this approach is its inability to explicitly treat 
three-particle states, such as 5s5p 2 . However, LCCS- 
DpT method produced the results for the 6pi/2 — 7s and 
6P1/2 ~ 7pi/2 Stark shifts of Tl, which is a similar 
Group IIIB system, in excellent agreement with the ex- 
periments [20 , [2~i| . For convenience, we will refer to the 
LCCSDpT approach as the CC method in text and tables 
below. 

The CI+all-order method (we also refer to it as CI+A11 
in text and tables below) was developed in Q. It allows 
us to treat In as a three-particle system, so all three elec- 
trons above the 4d 10 shell are considered valence. In this 
approach, the CC method is first used to accurately de- 
scribe core-core and core- valence correlation and to incor- 
porate them into the effective Hamiltonian. Therefore, 
the core-core and core-valence sectors of the correlation 
corrections for systems with few valence electrons will be 
treated with the same accuracy as in the all-order ap- 
proach for monovalent atoms. Then, the CI method is 
used to treat valence-valence correlations. Since the CI 
space includes only three electrons, it can be made essen- 
tially complete. This method allows to include correla- 
tion corrections to the wave functions in a more complete 
way than the CC approach. In particular it is capable to 
accurately account for the configuration mixing. 

The CI+all-order method yielded accurate wave func- 
tions for calculations of such atomic properties as life- 
times, polarizabilities, hyperfine structure constants, etc. 
for a number of divalent systems and Tl [si— Tl2l] . How- 
ever, the various types of the corrections to the effective 
dipole operator D c g are included in a more complete way 
in the CC approach at the present time. Therefore, both 
approaches are complementary and the difference in the 
results can serve as an estimate of the uncertainties. 

First, we compare the In energy levels calculated using 
the CC and CI+all-order methods with experiment (22j. 
The CC values include the part of the third-order energy 
not included by the LCCSD method and Breit interac- 
tion in second order. The CI+all-order method includes 
Breit interaction on the same footing as the Coulomb in- 
teraction, which would include some higher-order Breit 



TABLE I: Comparison of the CC and CI+all-order (labeled as 
"CI+ All" ) energies of In levels with experiment [23] • Three- 
electron binding energies are given in the first row. The en- 
ergies in other rows are given relative to the ground state. 
Corresponding relative differences of these two calculations 
with experiment are given in the corresponding columns la- 
beled "Diff." in %. The 5s 2 6p 2 P° /2 - 5s 2 6s 2 S 1/2 and 
5s 2 5d 2 £>5/2 — 5s 2 6p P3/2 transition energies are given in the 
last row. 



State 



5s 2 5p 



5s 2 6p 
5s 2 7p 

5s 2 6s 
5s 2 5d 

5s5p 2 
5s 2 7s 



Expt. CC Diff. CI+A11 Diff. 



2po 

"l/2 
2po 
r 3/2 
2 po 
-^1/2 
2po 

"s/2 
2po 
*l/2 
2po 
"3/2 

2 Sl/2 
D3/2 

2 D 5/2 

A/2 

-P3/2 
-P5/2 

Sl/2 



A(6pi /2 - 6s) 
A(5d 5/2 - 



P3/2J 



425060 

2213 2168 

31817 31468 

32115 31769 

38861 38513 

38973 38625 

24373 23862 

32892 32563 

32916 32754 
34978 
36021 
37452 

36302 35928 

7444 7606 

800 985 



2.02% 
1.10% 
1.08% 
0.90% 
0.89% 

2.10% 
1.00% 
0.49% 



1.03% 

-2.18% 
-23% 



425719 
2195 
31805 
32104 
38911 
39023 

24272 
32836 
32863 
35299 
36346 
37770 
36284 

7531 
759 



0.15% 
0.82% 
0.04% 
0.03% 
-0.13% 
-0.13% 

0.41% 
0.17% 
0.16% 
-0.92% 
-0.90% 
-0.85% 
0.05% 



16% 
16% 



corrections. The Breit correction is small; however, it sig- 
nificantly improves the accuracy of the 5s5p 2 4 Pj triplet 
splitting. Both calculations are carried out with l max — 6 
partial waves in all intermediate sums in the many-body 
expressions and include extrapolation for the contribu- 
tions of the partial waves with / > 6. 

In the CC calculations, extrapolation is carried out in 
second order using a separate code. This second-order 
calculation shows that the total contribution of the I > 6 
partial waves is close to the contribution of the single 
I = 6 partial wave. This empirical observation is used 
to estimate the contribution of the higher-order partial 
waves in the CI+all-order approach. We find that while 
both methods give energy levels in very good agreement 
with experiment, the CI+all-order results systematically 
agree better with the experimental values. 



III. POLARIZABILITIES 

Static polarizability of the state with total angular mo- 
mentum J, its projection M, and the energy E can be 
expressed as a sum over unperturbed intermediate states: 



a(J,M) = 2j2 



(J,M\D z \J n ,M)\' 

E n — E 



(1) 



where J n and E n are the total angular momenta and the 
energies of the intermediate states. 

The static polarizability a( J, M) can be conveniently 
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decomposed into scalar and tensor parts according to 

3M 2 - J(J+ 1) 



a = ao + ot2- 



J(2J-1) 



(2) 



where the scalar fan) and tensor (a 2 ) polarizabilities can 
be expressed as |23j ] 



a 



3(2J 



2 V l(./||^l|Jn)| 2 

1) ^ £L - E 



(3) 



and 



Q'2 



5J(2J- 1) 



1/2 



6(2J + 3)(2J + 1)(J + 1) 

For an open-shell atom, ao may be separated into a 
contribution from the valence electrons, ajj, contribution 
from the core electrons, a c , and a core modification due 
to the presence of the valence electrons, a vc . We calculate 
core and a vc terms using RPA in both approaches. The 
valence parts are calculated differently in the CC and 
CI+all-order methods. 

In the CC approach, the valence polarizability of the 
single-electron valence state |w) is calculated using the 
sum-over-states method: 



(k\\D\\ 



3(2j w + 1) 



k 



Ei. — E v 



(5) 



where j w is the total angular momentum of the state 
(k \\D\\ w) are the single-electron reduced electric-dipole 
matrix elements and sum over k runs over all intermedi- 
ate excited states allowed by the electric-dipole transition 
rules [13. 

In the CI+all-order approach, the valence part of the 
polarizability is determined by solving the inhomoge- 
neous equation of perturbation theory in the valence 
space, which is approximated as [25| 



{E-H eS )\^,M') = (D eS ) q \ 7 ,J,M), 



(6) 



where 7 incorporates all other quantum numbers except 
J and M. The wave function I7, M'), where M 1 = M+q, 
is composed of parts that have angular momenta of J' — 
J, J ± 1 . The construction of the effective Hamiltonian 
H e g using the all-order approach is described in [6j . The 
effective dipole operator D e g includes RPA corrections. 

While we do not use the sum-over-state approach in 
the calculations of the polarizabilities in the CI+all-order 
method, it is useful to calculate several dominant contri- 
butions to polarizabilities by combining the CI+all-order 
values of the El matrix elements and energies according 
to the sum-over-states formula ([5]) above. This allows us 
to compare dominant terms and total remainders calcu- 
lated by the CC and CI+all-order methods. 



TABLE II: Contributions to the 6s and 5pi/2 static polariz- 
abilities are given in a[] in columns labeled "ao" • The exper- 
imental energies (in cm -1 ) and the theoretical electric-dipole 
reduced matrix elements (in a.u.) used to calculate dominant 
contributions are listed in columns labeled "A_E" and "£>". 
The CC and CI+all-order electric-dipole matrix elements and 
the polarizability contributions are listed in columns labeled 
"CC" and "CI+AH", respectively. 



Contribution 


AE 




D 


a 






Expt. 


CC 


CI+A11 


CC 


CI+AU 


6s polarizability 


5pi/2 


-24:6(6 


1.911 


1.885 


-11.0(6) 


-10.7 


6pi/2 


7 AAA 


f\ 1 1 n 

O.l 1U 


O.Uol 


OO ( ) 




7f>l/2 


14488 


0.683 


0.648 


2.4(4) 


2.1 


8P1/2 


17454 


0.277 


0.265 


0.3(1) 


0.3 


(9 - 12)p 1/2 








0.17(3) 




(n > 12)p 1/2 








0.14(5) 




5P3/2 


-22160 


2.935 


2.899 


-28(1) 


-28 


6f»3/2 


7742 


8.529 


8.491 


687(24) 


681 


7P3/2 


14600 


1.131 


1.084 


6.4(1.0) 


5.9 


8P3/2 


17508 


0.489 


0.479 


1.0(2) 


1.0 


(9 - 12)p 3 /2 








0.6(1) 




(n > 12)pi /2 








0.6(2) 




Other 










23 


Core 








29.6(4.0) 


3.2 


VC 








0.0 


0.0 


Main 








1024.9 


1015.5 


Remainder 








31.1 


26.2 


Total 








1056(27) 


1042 


5pi/2 polarizability 


6s 


24373 


1.911 


1.885 


11.0(6) 


10.7 


7s 


36302 


0.548 


0.534 


0.61(3) 


0.57 


8s 


40637 


0.297 


0.297 


0.16(1) 


0.16 


(9- 12)s 








0.13(3) 




(n > 12)s 








0.7(2) 




5^3/2 


32892 


2.623 


2.577 


15.3(1.5) 


14.8 


6^3/2 


39049 


1.001 


0.865 


1.9(2) 


1.40 


7d 3/2 


41836 


0.537 


0.308 


0.5(1) 


0.17 


(8 - 12)d 3/2 








0.6(2) 




(n > 12)d 3/2 








6.1(4.0) 




Other 










31.6 


Core" 








29.6(2.5) 


3.22 


VC 








-5.0 


-0.15 


Main 








29.4 


27.7 


Remainder 








32.1 


34.6 


Total 








61.5(5.6) 


62.4 



"The uncertainty is the sum of the core and vc uncertainties. 



We compare the contributions to the 6s and 5pi/2 po- 
larizabilities ag calculated by two approaches in Table [TT1 
The CC and CI+all-order reduced electric-dipole matrix 
elements and the contributions to the polarizability ao 
are listed in columns labeled "CC" and "CI+AU" , respec- 
tively. The experimental energies listed in column AE 
are used to calculate the dominant contributions for con- 
sistency and to improve accuracy. We hnd generally very 
good agreement of the CC and CI+all-order results, with 
the exception of the 5pi/2 — 6^/2 and the 5pi/2 — 7(^3/2 
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TABLE III: Contributions to the 6pi/2 static polarizability 
are listed (in do) in columns labeled "ao" ■ The experimental 
energies (in cm" 1 ) and the theoretical electric-dipole reduced 
matrix elements (in a.u.) used to calculate dominant contri- 
butions are listed in columns labeled "A_B" and "D" . The CC 
and Cl+all-order matrix elements and the polarizability con- 
tributions are listed in columns labeled "CC" and "CI+AU", 
respectively. The contributions from the terms 6s, 7s, 8s, 
5^3/2, 6^3/2, and 7(4/2 are grouped together in row "Main". 



Contribution 


AE 




D 


«o 






Expt. 


CC 


CI+AU 


CC 


CI+AU 


6s 


-7444 


6.110 


6.081 


-367(12) 


-363 


7s 


4485 


6.289 


6.239 


645(5) 


635 


8s 


8820 


1.294 


1.317 


14 


14 


(9- 12) s 








5(1) 




(n > 12)s 








4(1) 




5d 3 /2 


1075 


10.095 


9.893 


6933(140) 


6659 


6^3/2 


7232 


6.470 


6.477 


423(65) 


424 


7d 3 /2 


10019 


2.861 


2.848 


60(7) 


59 


(8 - 12)d 3/2 








35(11) 




(n > 12)d 3/2 








35(6) 




Other 










81 


Core 








30(4) 


3.2 


Main 








7709(154) 


7429 


Remainder 








108(13) 


84 


Total 








7817(155) 


7513 


Recommended 








7817(300) 





cases. These matrix elements are small and have very 
large correlation corrections. The uncertainties in the 
CC contributions are evaluated using the well-defined ap- 
proach descried in detail in [5| . It involves calculation of 
the spread of four different CC calculations of increasing 
accuracy for each matrix element. The results labeled 
"Other" in the "CI+AU" column are obtained by sub- 
tracting the separately listed dominant terms from the 
final valence value. 

We note that the core contribution in the CC approach 
is substantially larger then the one in the Cl+all-order 
approach since the 5s shell is included in the core in the 
CC method, but not in the Cl+all-order one. The un- 
certainty of the CC core+vc term is evaluated as the 
difference of the DHF an RPA total core+vc values. The 
uncertainties in the tail are determined based on the dif- 
ference of the RPA and all-order values for terms with 
n = 9- 12. 

As an additional comparison between two approaches, 
we group eight contributions to the 6s polarizability 
(6s— (5— 8)^1/2.3/2) together and list them in row "Main". 
The contributions from six terms (6s, 7s, 8s, 5^3/2, 6e?3/2, 
and 7d 3 / 2 ) are grouped together in row "Main" for the 
5pi/2 polarizability. These main terms are subtracted 
from the totals to obtain the contributions from remain- 
ing terms. These terms are listed in rows labeled "Re- 
mainder" . The difference of the remainder part for the 6s 
polarizability appears to indicate lower value of the core 
polarizability when 5s 2 is included in the core. Overall, 



there is good agreement of both main and remainder part 
between the two calculations. The final values for the 6s 
and 5pi/ 2 polarizabilities, the 6s — 5j»i/ 2 Stark shift, and 
their uncertainties are discussed in the next section. 

Contributions to the 6pi/2 and 67)3/2 polarizabilities of 
In (in Oq) are listed in Tables IIIII and IIV1 respectively. 
The tensor polarizability a 2 of the 6p 3 / 2 state is given 
in Table IIVI These tables are structured in exactly the 
same way as Table HT1 The results of both CC and Cl+all- 
order calculations are given. The only difference is the 
listing of the 5s5p 2 4 P§/2 contribution to the 6p 3 / 2 po- 
larizability. The contributions of the 5s5p 2 A Pj state to 
all other polarizabilities considered here were found to be 
negligible. 

We find 3.9% difference between the CC and Cl+all- 
order values for both scalar 6pj polarizabilities and 3.1% 
difference between the values of the tensor a 2 polariz- 
ability for the 6^3/2 state. These differences are caused 
by the 2% difference in the values of the 6pi/ 2 — 5d 3 / 2 , 
6p 3 /2 — 5d 3 / 2 , and 6p 3 / 2 — 5d 5 / 2 matrix elements which 
dominate the 6pj polarizability values. This 2% dif- 
ference is consistent with the expected accuracy of the 
Cl+all-order method for these transitions. 

We evaluate the accuracy of the CC values to be on 
the order of 1%. The uncertainty evaluation is carried 
out differently for these three transitions owing to a con- 
vergence problem of the all-order equations for the bd 
states. We performed the calculations with three and 
five iterations in the LCCSD approximation and in the 
LCCSDpT approximation, and carried out the scaling 
procedure (described, for example, in Ref. (l[) using these 
four different starting points. The spread of the result- 
ing scaled values is 0.7%. Since the scaling estimates 
the dominant omitted corrections in such transitions, all 
other omitted corrections should not exceed this upper 
bound of 0.7%, resulting in total uncertainty estimate of 
1%. The CC and Cl+all-order values for all other con- 
tributions, including small remainders, were found to be 
in a good agreement. Because we expect the CC method 
to yield more accurate values of the 6s — 5d matrix el- 
ements, we take the CC values of the 6p polarizabilities 
as the final results. 

Most likely, the discrepancy with the Cl+all-order cal- 
culations is caused by an omission of the small corrections 
to the effective dipole operator in the Cl+all-order ap- 
proach. However, it might be possible that the 6s — hd 
matrix elements are affected by the small mixing of the 
even 5s 2 5d states with the 5s5p 2 4 Pj triplet, which is ac- 
counted for by the Cl+all-order method, but not the CC 
method. The weight (in probability) of the 5s5p 2 config- 
uration in the 5dj levels is 0.02 — 0.03, i.e., small but non 
negligible. Moreover, there is a very strong mixing of the 
nd configurations. Therefore, we take the difference of 
the Cl+all-order and CC results as the final uncertainty 
to account for the possible uncertainty due to such mix- 
ing. We note that 5s 2 6s state essentially does not mix 
with the 5s5p 2 configuration (its weight is only 0.0003), 
so 6s — 6p matrix elements are not affected by such mix- 
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TABLE IV: Contributions to the 6P3/2 scalar and tensor polarizabilities are listed (in Oq) in columns labeled "ao" and "0:2". 
The experimental energies (in cm -1 ) and the theoretical electric-dipole reduced matrix elements (in a.u.) used to calculate 
dominant contributions are listed in columns labeled "A_E" and "D" . The CC and Cl+all-order matrix elements and the scalar 
and tensor polarizability contributions are listed in columns labeled "CC" and "CI+ All" , respectively. The contributions from 



the terms 6s, 7s, 8s, 


5dj, 6dj, and 7dj are g 


;rouped together in 


row "Main". 








Contribution 


AS 




D 






C*2 






Expt. 


CC 


CI+A11 


CC 


CI+A11 


CC 


CI+A11 


fi« 


-7742 


8.529 






-Oil 




on 


7s 


4187 


9.413 


9.338 


774(7) 

1 1 Try 1 y 


762 


-774(7) 


-762 


8s 


8522 


1.797 


1.830 


14 


14 


-14 


-14 


(9 - 12)s 








5 




-5 




(n > 12)s 








4(1) 




-4(1) 




oa 3/2 


777 


4.511 


4.420 








7^fi 


64/2 


6933 


3.157 


3.124 


53(7) 


51 


42(6) 


41 


1 u '3/2 


9721 


1.350 


1.309 


7(1) 


6 


5(1) 


5 


(8 - I2W3/0 








4(1) 




3(1) 




(n > 12)^/2 








4(1) 




3(1) 




5d 5 /2 


800 


13.577 


13.296 


8426(170) 


8080 


-1685(34) 


-1616 


6d 5 /2 


uyoo 


Q 91 8 


9.031 


445(58) 


427 


-89(12) 


-85 


7d 5 /2 


9747 


4.003 


3.930 


60(7) 


58 


-12(1) 


-12 


(8 - 12)d 5/2 








35(7) 




-7(1) 




(n > 12)d 5/2 








33(9) 




-7(1) 




5s5p 2 4 P 5 /2 


5337 




2.337 




37 




-7 


Other 










79 




-13 


Core 








30(4) 


3 






Main 








10393(180) 


9979 


-1416(42) 


-1367 


Remainder 








114(12) 


119 


-16(2) 


-20 


Total 








10506(180) 


10098 


-1432(42) 


-1387 


Recommended 








10500(400) 




-1432(45) 





ing. A high-precision measurement of the 6^3/2— 6s Stark 
shift should resolve this question. The final values for 
the 6^ 3/2 polarizabilities are listed in the last rows of 
Tables |ml and El 



IV. DISCUSSION OF THE RESULTS AND 
THEIR UNCERTAINTIES 

While we have estimated the uncertainties of the CC 
calculation, it is possible to improve our evaluation of 
the uncertainties by comparing the CC and Cl+all-order 
results. As we have described above, these two high- 
precision approaches include somewhat different higher- 
order effects. The Cl+all-order calculations include the 
valence- valence correlation corrections to the wave func- 
tions very precisely, as indicated by excellent agreement 
of the respective energies with experiment. On the 
other hand, the effective dipole operator D e g includes 
only RPA corrections in the Cl+all-order method at the 
present time, omitting the structure radiation, normal- 
ization, and other small corrections. These corrections 
to the electric-dipole matrix elements are included in the 
CC method. 

In the polarizability calculations, we use the experi- 
mental energies for the main terms. Therefore, more ac- 
curate method of calculating individual matrix elements 



is somewhat more important for the 6s polarizability. In 
the framework of the CC method, the 5s shell is included 
in the core and we have the large core polarizability. It 
leads to increasing the total uncertainty of the polariz- 
ability because the core polarizability is calculated with 
less accuracy But it cancels out when the Stark shift of a 
transition is calculated. The Cl+all-order method treats 
contributions with high n with better accuracy, since it 
is done by solving the inhomogeneous equation instead of 
using the sum-over-states method. This is not significant 
for the 6s polarizability where such tail contributions are 
small, but is important for the 5p±/2 polarizability. 

In summary, the CC and Cl+all-order methods to- 
gether include all correlation corrections that are ex- 
pected to be dominant for the present polarizability cal- 
culations. Therefore, the uncertainty can be approxi- 
mated as the difference 6a = |a(CC) — cv(CI+all)| of the 
CC and Cl+all-order results. All other omitted higher- 
order corrections are expected to be smaller than already 
included ones, therefore, we can also assume that they do 
not exceed 6a. Therefore, we take our final uncertainty 
in the polarizabilities and their difference to be \/25a, 
calculated separately for each of the properties. We list 
our final values of the 6s and 5pi/2 polarizabilities and 
their difference Aao(6s — 5pi/ 2 ) in Table IVl in a.u. 

We take the CC results to be the final values since CC 
method accounts for more higher-order corrections to the 
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TABLE V: Final values of the 6s and 5pi/2 polarizabilities and their difference Aao(a.u.). Determination of the reduced 
electric-dipole 6s — 6pj matrix elements (in a.u.) and 6pj lifetimes (in ns) from the combination of the recently measured Stark 
shift [l| and theoretical values. The quantity C is the value of Aao(6s — 5pi/2) with the contribution of 6s — 6pj transitions 
subtracted out. The results are compared with other theory and experiment. 





a (6s) 


ao(5pi/ a ) 


Aq (6s - 5pi/ 2 ) 


C 


D(6s-6pi /2 ) D(6s - 6p 3/2 ) 


T(6pi/a) 


T(6p 3 /2) 


cc 


1056 


61.5 


995 


-59.6 


6.126 8.551 






CI+A11 


1042 


62.4 


980 


-65.2 


6.141 8.575 






Final 


1056(20) 


61.5(1.3) 


995(21) 


-59.6(7.8) 


6.126(24) 8.551(34) 


63.77(50) 


58.17(45) 


Expt. [13] 






1000.2(2.7) 










Expt. [28] 




69(8) 












Theory [27] 




61.48 












Theory [26] 




62.0(1.9) 












Theory [29] 












63.8(8) 


58(1) 



El matrix elements that dominate the 6s polarizability. 
For consistency, CC value of the 5pi/ 2 polarizability is 
used when calculating the final value for the 6s — 5pi/ 2 
Stark shift. Our final result is in excellent agreement 
with recent high-precision measurement of the 6s — 5pi/ 2 
Stark shift [13J, which allows for benchmark compari- 
son of the theoretical approachers. We find that the CC 
value is closer to the experimental measurement than the 
CI+all-order result. Our calculated polarizability of the 
5pi/2 state is in excellent agreement with recent coupled- 
cluster single, double, and perturbative tripes excitation 
[CCSD(T)] calculations of Refs. [H, S3- We note that 
our implementation of the coupled-cluster method dif- 
fers significantly from that of Refs. (2|| H3] (see recent 
review [24| for details). 



V. DETERMINATION OF THE LIFETIMES 

Recent precision measurement of the 6s — 5pi / 2 Stark 
shift [l3| can be combined with the present calculations 
to obtain very accurate lifetimes of the 6pi/ 2 and 6p3/ 2 
states. This is possible since the 6s — 6pj matrix elements 
overwhelmingly dominate the values of this Stark shift as 
illustrated by Table HT1 Separating the 6s — 6pj contribu- 
tions (see Eq.©), we write the Aa (6s — 5pi/ 2 ) Stark 
shift as 



Aa (6s-5p 1/2 ) =BS + C, 



where 



1/1 R 2 

~ 3 \E{6 Pl/2 ) - E(6a) E{6p 3/2 ) - E{6s) 



(7) 



(8) 



S = D 2 is the 6s — 6p!/ 2 line strength, R is the ratio of 
the D(6s — 6p3/ 2 ) and D(6s — 6pi/ 2 ) reduced El matrix 
elements, and the term C contains all other contribu- 
tions to the Stark shift. We calculate the ratio R to be 
1.396(1). The uncertainty (0.001) is very small since the 
ratio R is very insensitive to different corrections. Using 
the results presented in Table HTl we see that the CC and 
CI+all-order methods give R equal to each other up to 



fourth significant figure. Substituting the ratio R and 
the corresponding experimental energies in Eq. (|5J| gives 
B = 28.24(3) a.u. 

The values of C in the CC and CI+all-order methods 
are obtained from the results presented in Table [TT] and 
are equal to -59.6 a.u. and -65.2 a.u., respectively (see Ta- 
ble|V|. Adding the relevant uncertainties from CC calcu- 
lations of Table [II] in quadrature leads to the uncertainty 
SC = 4.9 a.u. This number is consistent with the differ- 
ence of the CC and CI+all-order values, 5.6 a.u. We note 
that the uncertainty of the core term does not contribute 
to SC, since this term is canceled when the Stark shift is 
calculated. The 5pi/ 2 a vc term does contribute, and its 
uncertainty is 1.3 a.u. To evaluate the final uncertainty 
in C we use the same rule as for the polarizability: mul- 
tiply the difference of the CC and CI+all-order values by 
V2 to account for other small uncertainties not included 
in our consideration. Again, we assume that they cannot 
exceed the difference of the CC and CI+all-order values. 
The final value for C is presented in Table |Vj 

There are three sources of the uncertainties contribut- 
ing to the uncertainty in D(6s — 6p!/ 2 ): uncertainties in 
C, B, and experimental values of Aao- For convenience, 
we calculate first the uncertainty in line strength S using 
formula: 



SS 



^/(SC) 2 + (SAa ) 2 + (SSB) 2 



(9) 



The relative uncertainty in I? is a half of the relative 
uncertainty in S. The uncertainty in B is negligible. The 
final values of the matrix elements are listed in Table [V] 
The lifetimes of the 6pi/ 2 and 6p3/ 2 states are obtained 
using the formula r a = 1/A a b since there is only one El 
decay channel for each state. The transition rate A a b is 
given by 



A„ 



2.02613 x 10 



18 



S a b 



XL 



2J a + l 



(10) 



where the transition wavelength A a {, is in A . The relative 
uncertainties in the lifetimes are twice that of the relative 
uncertainties of the corresponding El matrix elements. 
The final values are given in Table [V] 
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To simplify the extraction of the 5dj lifetimes from fu- 
ture experimental value of the 6p 3 / 2 — 6s Stark shift, we 
evaluated the sum of all contributions to the 6pi/2 and 
6p 3 / 2 polarizabilities with the 5d — 6p terms excluded. 

These quantities, which we designate as C(6pj), are ob- 
tained from the data in Tables IIIII and IIVI We note that 
C(6pi/ 2 ) and C(6p 3 / 2 ) refer to the contributions to the 
polarizabilities, rather than their differences as in Eq. ([7]). 

We find that the CC and CI+all-order results are 
very close together, and are well within the uncertainty 
estimates of the CC data. The CC values (in a.u.) 
are C(6p 1/2 ) = 884(68), C (6p 3 / 2 ) = 1123(60), and 

C2(6p 3 / 2 ) = —514(19), where So and 5*2 are related to the 
scalar and tensor polarizabilities, respectively. The corre- 
sponding CI+all-order values are 854, 1098 and —506 (in 
a.u.). The ratio of the 6pi/2 — 5d 3 /2 and 6p 3 / 2 — 5<i 3 / 2 ma- 
trix elements is 2.238(4) and the ratio of the 6p 3 / 2 — 5d 5 / 2 
and 6p 3 / 2 — 5<i 3 /2 matrix elements is 3.0095(16). 

VI. CONCLUSION 

We carried out a first systematic comparison of the 
linearized coupled-cluster and CI+all-order method us- 
ing the polarizabilities of the low-lying states of the In 
atom as a benchmark testing case. We find that the 
CI+all-order method produces more accurate data for 
the energy levels. It appears that the CC data for the El 
matrix elements are somewhat more accurate due to more 
complete inclusion of the small higher-order corrections 



to the matrix elements in the cases where relevant con- 
figuration mixing of trivalent states is negligible. This is 
an additional motivation to incorporate such corrections 
into the CI+all-order formalism at the all-order level in 
the future. 

Our result for the 6s — 5pi/2 Stark shift is in excellent 
agreement with the recent high-precision experiment [l3[ . 
We also provide predictions for the polarizabilities of the 
Qpi/2 and 6p 3 / 2 states. A precise experimental measure- 
ment of the 6p 3 /2 — 6s Stark shift proposed in 13 would 
be a good test of our calculations. This will also provide 
an excellent test of the theoretical approaches. Combin- 
ing the present calculations with the experimental Stark 
shift data allows very accurate extraction for the lifetimes 
of the low-lying In states. 
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